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Abstract 

The generalized non-extensive statistics proposed by Tsallis have been successfully utilized in 
many systems where long range interactions are present. For high density quark matter an attrac- 
tive long range interaction arising from single gluon exchange suggests the formation of a diquark 
condensate. We study the effects on this color superconducting phase for two quark flavors due 
to a change to Tsallis statistics. By numerically solving the gap equation we obtain a generaliza- 
tion of the universality condition, ~ 3.52 and determine the temperature dependence of the 
gap. For the Tsallis parameter g ~ 1 the specific heat is exponential becoming more linear as q 
increases. This suggests that for larger values of q s-wave color superconductors behave like high 
Tc superconductors rather than weak superconductors. 
PACS: 12.38.Mh,21.65.QR 
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INTRODUCTION 



Nonextensive thermostatistics 



QflflQ 



have been utilized with success in connection 



with a number of problems both in the classical [5, 6, ^jls!,!^, 11 1 and the quantum regimes 
3, Q, Q, H, Il6l . llTI ]. They are thought to be relevant for the study (among others) of: 
systems described by non linear Fokker-Planck equations [5]; systems with a scale-invariant 



occupancy of phase space Gj; non equi 



ibrium scenarios involving temperature fluctuations 



and systems with interactions of long range 



systems exhibiting weak chaos 
relative to the system's size 11 1. 

n 

Recently Hagedorn[18| statistical theory of the momentum spectra produced in heavy ion 
collisions has been generalized using Tsallis statistics to provide a good description of e^e~ 
annihilation experiments [20|. Furthermore, Walton and Rafelski[2l| studied a Fokker- 
Planck equation describing charmed quarks in a thermal quark-gluon plasma and showed 
that Tsallis statistics were relevant. These results suggest that perhaps Boltzmann-Gibbs 
statistics may not necessarily be adequate in the quark-gluon phase. 

The existence in QCD of single gluon exchange betwee n q uarks of different color generates 
an attractive interaction in the color-antitriplet channel 22|] which leads unavoidably to the 



occurrence of a particle-particle condensate 



23 



241 ] or color superconductivity in dense quark 



matter at low temperatures. Color superconductivity differs from normal superconductivity 
which is described by BCS-like theories which usually assumes an interaction which is heavily 
screened. For a point-like four fermion coupling there is no correlation between the initial and 
outgoing momenta of the interacting fermions. Hence the BCS gap is constant with respect 
to the momentum as long as the momenta are near the Fermi surface. In QCD, however, 
scattering of quarks through single gluon exchange exhibits a logarithmic divergence for 
forward scattering. As a result the gap is not exponential in ^, where g is the coupling 
constant for QCD, as in BCS but only in - 
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271, 



281]. Hence the gap is no longer 



constant as function of the momentum, even around the Fermi surface. The logarithmic 
divergence arises because the static magnetic interactions are not screened in dense quark 
matter at low temperatures. 

The absence of a fully screened interaction in color superconductivity suggests that per- 
haps it might be useful to examine the role that the choice of statistics plays in this phe- 
nomenon. We therefore have solved the gap equations using the generalized thermostatistics 
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of Tsallispl] rather than that of Boltzmann Gibbs. In the hmit where TsaUis parameter goes 
to one we recover the Boltzmann Gibbs results. We show that the universality condition, 
which characterizes BCS-like superconductivity, is easily generalized if TsaUis statistics are 
used. Furthermore, calculation of the specific heat as a function of the Tsallis parameter 
suggest that as the parameter increases s-wave color superconductors may behave more like 
high Tc superconductors than weak superconductors. 



FORMALISM 



Rischke and Pisarski (RP)have developed a formalism for the study of color supercon- 
ductors in the weak coupling regime. This formalism is valid for systems involving large 



quark densities and low temperatures 
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29|. One advantage of this formalism is that the 



temperature dependence of the gap equation is completely transparent. As we will discuss, 
it is the temperature dependence of the gap equation that is modified when one introduces 
Tsallis statistics. 

In general two fundamental color triplet states can pair up into a symmetric color sextet 
and an antisymmetric color antitriplet. The antitriplet channel of quark scattering is attrac- 
tive and leads to the color condensate. It is well known that these condensates can not form 
a color singlet. Therefore, SU{3)c is broken in this phase. From a single quark flavor one 
cannot form a color antitriplet. The simplest possibility then, and the one we will focus on, 
is the case of two quark flavors forming a J=0 condensate. This is the so-called 2SC phase. 

In the RP formalism the gap equation is derived in momentum space using the Hard 
Dense Loop (HDL) approximation. In the HDL approximation the gluon propagator re- 
ceives a mass m ^ gfi where g is the QCD coupling and fi is the quark chemical potential. 
This approximation generates a hierarchy of scales fi ^ nig ^ (po where 0o is the gap at 
temperature T=0. In this approximation the gap equation is 



tanh(-^)^ln(— 1^)0, (0.1) 



where b = 2567r^(p^)^/2, 5 is a cutoff scale, and = ^/{q — yu)^+ | p is the quasi- 
particle excitation energy. This integral equation can be turned into a soft integral by the 
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approximation 

I M l f^^i ^ H-)0{q -k) + H^)e{q - k). (0.2) 
^ I '^g ^fc I <^fc 

Note that the resuhing energy gap is a function of momentum k. This is different from 
BCS superconductors and is due to the long range nature of the color magnetic interaction. 
Following RP we introduce g = and change variables to y = ln( ^_^^^^ ), x = ln( ^_^^^^^ ) 
and X* = In(^). This yields the following form for the soft integral equation: 

m ^ f[x r dytanhC-^my) + f dyy tanhC-^ (^)]. (0.3) 



2T J\n{bn/5 2T 

For T=0 this equation can be solved algebraically to give 

TT 



(x) = 26/i exp(— — ) sin(5fx) (0.4) 
^9 



at lowest order in g 



At nonzero temperature the gap equation must be solved numerically. Numerical solu- 
tions of Eqs. flO.ll) and (10.31) for the 2SC model have been studied extensively by Zakout et. 



al. 30|. Zakout et. al. obtained numerical solutions for the full gap equation (for a large 
coupling, g > 2.0) as well as the approximate (Eq. flO.ll) ) and soft integral (Eq. fl0.3l) ) forms 
for the gap for small values of g. Comparing to the the full gap equation, it was shown that 
the approximations were valid to leading order. The gap equation has a sharp peak at the 
Fermi surface, i.e. at A; = /i in the weak coupling regime. Zakout et. al. found that it was 
essential to use the soft integral form of the gap, Eq. (10.31) . in order to properly secure the 
singularity. For the gap equation obtained using Tsallis statistics, our numerical calculations 
lead us to the same conclusion. 

To determine the critical temperature, Tc, one assumes (f){x,T) ^ (/)(T)^^^ where 0o is 
the gap at T=0. With this assumption the gap equation at the Fermi surface becomes 

l^fT ,,,ta„h(iM)^. (0.5) 



RS have shown that this yields the standard BCS universality result 23| 



^ ^ 3.52. (0.6) 

One important issue that we address in the following section is whether universality can be 
maintained using Tsallis statistics. 



COLOR SUPERCONDUCTIVITY WITH TSALLIS STATISTICS 



We now investigate the effect on the gap equation from introducing generahzed Tsalhs 
statistics. The starting point is the non-extensive entropy proposed by Tsalhs: 

S, = fcS^^^Mi (0.7) 

where w is the total number of microstates, p is the associated probabilities, and q is a real 
number. By taking q ^ 1 one recovers the usual Boltzmann-Gibbs (BG) expression for 
entropy. This results in a modification of the usual Fermi-Dirac distribution. Working in 
units where k=c=h=l the generalized Fermi distribution becomes 

f, = 1 • (0.8) 

[l + (g-l)e/T]^ + l 

This results in a modification of the temperature dependence of Eq. (10.31) . In Eq. (10.31) the 
usual BG expression for the distribution function has already been assumed. Replacing 
the usual Fermi distribution with Eq. (10.81) leads to the following modification in the gap 
equation: 

2r (l + l/r(5-l)e,)— +1 

As discussed earlier it is essential to use the soft integral approximation when numerically 
solving the gap equation in the weak coupling regime. The generalized gap equation is 

i^ff ,,,(i + Vr(,-iK)A-i^(,.o)^ ^^^^.^ 

Jin{bt,/5) (l + l/T{q-l)eq)~ + 1 (Po 
For g we use the running coupling relation 

9 87r2 1 ^ ^ 

where I3q = [IINq — 2Nf)/3, Nq and Nf are the numbers of quark colors and flavors 
respectively, and A = 400MeV. For the integration cutoff we take 6 ~ rrig where 



< = Nf^-X + (Nc + ^)^ (0.12) 
In the regime we are considering the temperature dependence of Eq. (I0.12p may be ignored. 
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Clearly, the universality condition in its canonical form will not hold since Tc is now a 
function of q. However, as has already been noted in the case of high Tn BCS superconduc- 
tors, there is a reasonable generalization of the universality condition 141 ]: 

^^3.52 (0.13) 
qTc ^ ' 

Setting 0(g) = in Eg. (10. 101) and taking g ^ 0.39 we numerically solve the gap equation 

for Tc for different values of q. FIGH] shows both 2(j)o/Tc and 2(f)Q/qTc as a function of 

q. For reasonable values of q the generalized universality condition shows no appreciable 

deviation from 3.52. It is also possible that a better fit could be obtained by replacing q 

with some function of q. 



9 J ^ ^ ^ ^ ^ ^ 
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FIG. 1: 2(Pq/Tc (solid) and 2(/)o/(?Tc (dotted) as a function of q 

We next consider the temperature dependence of the gap equation Eq. OO.lOp . As men- 
tioned earlier, the gap is a function of the quasiparticle momentum k. In order to get 
quantitative results we evaluate the gap equation at the Fermi surface, k = ^. Although q 
can in principle be any real number, applications of Tsallis statistics typically suggest a q 



2^o/qT, 
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value not significantly greater than 1. Fig. [2] shows the temperature dependence of the gap 
with Tsallis statistics. In each case q is determined by fixing the ratio 20o/^c- 




FIG. 2: 2(t){T)/Tc (solid) and 2(t){T)/qTc (dashed) as a function of q 

For q=1.178, 0(r) closely resembles the usual behavior of the gap, while increasing 
q results in a flatter curve. By measuring the temperature dependence of the gap for 
superconducting quark matter one would clearly be able to differentiate the cases of Tsallis 
and Boltzmann-Gibbs statistics. However, it is likely easier to measure the specific heat of 
the quark matter. This can be calculated using the usual thermodynamic relation 

C^T% (0.14) 

where the entropy is given by 

S, = ^[ fd{q - ^)(/(6) - f{ey) + fd{q - /i)((l - f{e)) - (1 - f{e)Y)] (0.15) 
Q ^ ^ Jo Jo 

FIG. |3] shows the resulting specific heats for the same values of q considered in FIG. El 
For g 1, C displays a more exponential like behavior which is indicative of a weak 



superconductor 



23 



3l| . A more linear dependence on temperature occurs for larger values 
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FIG. 3: C /Cmax where Cmax is the maximum value of the specific heat 

of q suggesting that for these values of the TsaUis parameter s-wave color superconductors 
behave more like high Tc superconductors jl4 ] . 

A generalization of s-wave color superconductivity to include non-Boltzmann Gibbs en- 
tropic measures has been given. Introducing Tsallis statistics provides a simple means of 
studying the behavior of color superconductors in the limit of weak and high temperature 
BCS superconductors. Signatures are most easily identified by looking at the specific heat. 
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